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Abstract: The U, V method of estimation provides unbiased estimators or 
predictors of random quantities. The method was introduced by Robbins Q 
and subsequently studied in a series of papers by Robbins and Zhang. (See 
Zhang Practical applications of the method are featured in these papers. 
We demonstrate that for one U function (one for which there is an important 
application) the V estimator is inadmissible for a wide class of loss functions. 
For another important U function the V estimator is admissible for the squared 
error loss function. 

1. Introduction 

The U, V method of estimation was introduced by Robbins Q • The method applies 
to estimating random quantities in an unbiased way, where unbiasedncss is defined 
as follows: The expected value of the estimator equals the expected value of the 
random quantity to be estimated. More specifically, suppose Xj, j = 1, . . . ,n, are 
random variables whose density (or mass) function is denoted by fxii. x i\^i)- I n this 
paper we consider estimands of the form 

n 

(1.1) s(x,e) = Y / u*(x J ,e J ), 

where X = (X\, . . . , X n )' and 6 = (9%, . . . , 8 n )' . An estimator, V(X) is an unbiased 
estimator of S if 

(1.2) E e V(X) = E e (S(X,0)). 

Of particular interest in applications are estimands of the form U* (Xj ,6j) = 
U(Xj)6j, where U(-) is an indicator function. Robbins [|| offers a number of ex- 
amples of unbiased estimators using the U, V method. Zhang [|| studies the U, V 
method for estimating S and provides conditions under which the "U, V" estima- 
tors are asymptotically efficient. Zhang 0] then presents a Poisson example that 
deals with a practical problem involving motor vehicle accidents. 

In this note we demonstrate that for many practical applications the U, V es- 
timators are inadmissible for many sensible loss functions. In particular, for the 
Poisson example given in Zhang 0] , for the U function given, the V estimator is in- 
admissible for any reasonable loss function, since the estimator is positive for some 
X when S = no matter which 9 is true. 

Previously, Sackrowitz and Samuel-Cahn [1] showed that the U, V estimator of 
the selected mean of two independent negative exponential distributions is inad- 
missible for squared error loss. 
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In the next section we examine examples in which S functions based on simple U 
functions are estimated by inadmissible V functions. For other simple U functions 
the resulting V estimators are admissible for squared error loss. These later results 
will be presented in Section 3. 



2. Inadmissibility results 



Let Xj, j = 1, . . . , n, be independent random variables with density f Xi (xi\9i). Let 
U*(Xj,9j) = U{Xj)dj, where, for some fixed A>0, 



(2.1) U(Xj) 



l,ifXj < A, 



0, if Xj > A. 
Consider the following four distributions for Xj. 

(2.2) Poisson f x {x\6) =e-°6 x /x\ {9 > 0,x = 0, 1, . . .), 

(2.3) Geometric f x {x\9) = (1 - 6)6 X (0 < 9 < l,x = 0, 1, . . .), 

(2.4) Exponential f X {x\9) = (l/9)e~ x/e (9>0,x>0), 

(2.5) Uniform Scale f X (x\6) = 1/9 {Q<x<9.9> 0). 

Let W(t), t > be a function with the property that W(0) = and W(t) > for 
t > 0. Consider loss functions 



(2.6) W(a,S) = W{a- S), 

for action a. 

For the distributions in (2.2), (2.3), (2.4), (2.5), Robbins 0] finds unique unbiased 
estimators V(Xj) for U{Xj)6j. 

Theorem 2.1. Let Xj, j = 1, . . . , n, be independent random variables whose distri- 
bution is (2.2) or (2.3) or (2.4) or (2.5). Consider the loss function given in (2.6). 
Let U(Xj) be as in (2.1). Then the unbiased estimator V (X) = Y)j—\ V(Xj), where 
V{Xj) is the unbiased estimator of U{Xj)9j , is inadmissible for S given in (1.1). 

Proof. The idea of the proof is easily seen if n = 1. However for n > 1 it is 
instructive to see how much improvement can be made. The proof for n — 1 goes 
as follows: Let X\ be X and 9\ be 9. The V(X) estimators for the four cases are 
given in Robbins @|. For the Poisson case V(X) = U(X - 1)X (V(0) = 0). Now let 
[A] denote the largest integer in A less that A. Then V([A] + 1) = [A] + 1, whereas 
S = U([A] + 1)0 = 0. 
If 

_ / V(X), all X except X = [A] + 1, 
V {X) -\0, X=[A] + 1, 

then clearly V*(X) is better than V(X) since W(V*([A] + 1) - S) = for V* 
and W(([^4] + 1) — S) > for V. For the case of arbitrary n, S = whenever all 
Xj > ([A} + 1) whereas V(X) ^ whenever at least one Xj = ([A] + 1). If all 
Xj = ([A] + 1), then V = n([A] + 1). Clearly if V* = at such X, V* is better 
than V. 

For the geometric distribution when n = 1, V(X) = Ylif=^U{i) (V(G) = 0). 
Note S = for X > [A] + 1 but V = [A] + 1 for all such X. Again if V* = V for 
X < [A] and V* = for X > [A] + 1, V* is better than V. The case of arbitrary 
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Table 1 

Improvement in risk for squared error loss function 



n 



A 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


1 


1.083 


1.872 


2.190 


2.148 


1.902 


1.575 


1.243 


0.947 


0.701 


0.508 


3 


3.126 


4.763 


5.086 


4.626 


3.831 


2.982 


2.220 


1.599 


1.122 


0.771 


5 


5.782 


8.268 


8.419 


7.364 


5.894 


4.447 


3.216 


2.253 


1.539 


1.031 


7 


8.934 


12.268 


12.113 


10.328 


8.083 


5.976 


4.242 


2.919 


1.961 


1.292 


9 


12.511 


16.694 


16.120 


13.490 


10.388 


7.568 


5.299 


3.600 


2.389 


1.556 



n is even more dramatic than is the Poisson case with S — if all Xj > [A] + 1 
whereas V ^ on such points. 

For the exponential distribution when n = 1, V(X) = L U(t)dt = X if X < A, 
and V(X) = A if X > A. For arbitrary n, S = whenever all Xj > A, whereas 
V{X) ^ on such points. 

For the scale parameter of a uniform distribution with n = 1, V(X) — XU(X) + 
J* U(t)dt which becomes 2X if X < A and A if X > A. Hence as in the previous 
case, for arbitrary n, S = whenever all Xj > A whereas V(X) ^ on such points. 
This completes the proof of the theorem. □ 

Remark 2.1. Theorem 2.1 applies to the Poisson example in Zhang 

Remark 2.2. // the loss function in (2.6) is squared error then the amount of 
improvement in risk ofV* over V depends on n, A, and 9. It can be easily calculated. 
For the case where all the components of 9 are equal and each 6i, i = 1, . . . ,n is 
set equal to [A] + 1 the amount of improvement is equal to 

e:Li (m + i)) 2 c? e -^)({A] + 

Table 1 offers the amount of improvement for n = 1(1)10 and for values of 
A = 1,3,5,7,9. We observe as n gets large the amount of improvement becomes 
smaller. Also for small n as A gets large, improvement gets large. Such observations 
are consistent with the asymptotic efficiency of the U, V estimator as n — > oo and 
with Sterling's formula. 

Remark 2.3. Theorem 2.1 also holds for predicting 

n 

s * = J2 Y * u (Xj), 

3=1 

where Yj has the same distribution of Xj but is unobserved. 




3. Admissibility results 

In this section we consider the case 

(3.i) vto) = {°i£%ii "4-0; j = i) 



Also we consider a squared error loss function. 
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Theorem 3.1. Suppose Xj are independent with Poisson distributions with pa- 
rameter Xj. Then V(X) is an admissible estimator of S(X,\) for squared error 
loss. 



Proof. Let n = 1 and recall V(Xi) = U(X 1 - 1)X U V(0) = 0. Then 
V(X) 



0, forXi = 0,1,...,L4] + 1, 
X u torXi > [A] + l, 



while 

U*{X 1 ,X 1 ) = U(X 1 )X 1 



0, X X <[A], 
Xx,X x >[A] + 1 



Since U*(Xi, Ai) = for X\ < [A], any admissible estimator of U* (Xi, Ai) must 
estimate for X\ < [A] as V(Xi) does. □ 

At this point we can restrict the class of estimators to all those which estimate 
by the value for all X ± < [A]. For [Xi] > [A] + 1, U*(Xi, Ai) = Ai and we have a 
traditional problem of estimating a parameter Ai. Now we can refer to the proof of 
Lemma 5.2 of Brown and Farrell [H to conclude that any estimator that can beat 
V(X) would have to estimate at X\ — [A] + 1. Furthermore for the conditional 
problem given X% > [A] + 1, it follows by results in Johnstone 0] that X% is an 
admissible estimator of Ai. 

For arbitrary n the proof is more detailed. We give the details for n — 2. The 
extension for arbitrary n will follow the steps for n — 2 and employ induction. For 
n = 2, suppose V(Xi) + V(X<2.) is inadmissible. Then there exists S*(X\,X2) such 
that 



oc oc 



E E (^(^i) + ^2)-C/(a ; i)Ai-{/( a ; 2 )A 2 )'A^A^ e - Al - A V a ; 1 



\x 2 \ 



X\— X2— 



(3.2) > {5*{x 1 ,x 2 )-U{x 1 )X 1 -U{x 2 )X 2 ) 2 X x 1 1 X x 2 i e-^-^/x 1 \x2\ 



xi= xi— 



for all Ai > 0, A 2 > 0, with strict inequality for some Ai and X 2 . 

Now let A2 — * 0. Then by continuity of the risk function, (3.2) leads to 

(3.3) £<jV(*i) - ^i)Ai) 2 | > E^(5*(X U 0) - C/(X 1 )A 1 ) 2 |. 

Since V(Xi) is admissible for U(Xi)Xi, the case n = 1, (3.3) implies that V(Xl) = 
<5*(Xi,0). At this point we do as in Brown and Farrell [H by dividing both sides 
of (3.2) by A 2 . Reconsider (3.2) but now we can let the sum on x 2 run from 1 
to 00 since V{X\) = S*(Xi,0). Again let A 2 — > and this leads to V{X\) = 
6*(X±, 1). Repeat the process for X 2 = 0, 1, . . . , [A] + 1. Furthermore by symmetry 
V(X 2 ) = S*(0,X 2 ) = ••• = S*([A] + l,X 2 ). Thus V(X\) + V(X 2 ) = 5*(Xi,X 2 ) 
on all sample points except the set B = (X\ > [A] + 2,X 2 > [A] + 2). Here 
V(X±) + V(X 2 ) = Xi + X 2 and S = Ai + X 2 . We consider the conditional problem 
of estimating Ai + A2 by X% + X 2 given X 6 B. Clearly when Ai = X 2 = X no 
estimator can match, much less beat the risk of A1+X2 for this conditional problem 
since X\ + X 2 is a sufficient statistic, the loss is squared error, and X\ + X 2 is an 
admissible estimator of 2A. Thus S*(Xi,X 2 ) — V(Xi) + V(X 2 ) on the entire sample 
space proving the theorem. 
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